Correlated substitutions reveal SARS-like coronaviruses recombine frequently with a diverse set of structured gene pools.

Quantifying SARS-like coronavirus (SL-CoV) evolution is critical to understanding the origins of SARS-CoV-2 and the molecular processes that could underlie future epidemic viruses. While genomic analyses suggest recombination was a factor in the emergence of SARS-CoV-2, few studies have quantified recombination rates among SL-CoVs. Here, we infer recombination rates of SL-CoVs from correlated substitutions in sequencing data using a coalescent model with recombination. Our computationally-efficient, non-phylogenetic method infers recombination parameters of both sampled sequences and the unsampled gene pools with which they recombine. We apply this approach to infer recombination parameters for a range of positive-sense RNA viruses. We then analyze a set of 191 SL-CoV sequences (including SARS-CoV-2) and find that ORF1ab and S genes frequently undergo recombination. We identify which SL-CoV sequence clusters have recombined with shared gene pools, and show that these pools have distinct structures and high recombination rates, with multiple recombination events occurring per synonymous substitution. We find that individual genes have recombined with different viral reservoirs. By decoupling contributions from mutation and recombination, we recover the phylogeny of non-recombined portions for many of these SL-CoVs, including the position of SARS-CoV-2 in this clonal phylogeny. Lastly, by analyzing >400,000 SARS-CoV-2 whole genome sequences, we show current diversity levels are insufficient to infer the within-population recombination rate of the virus since the pandemic began. Our work offers new methods for inferring recombination rates in RNA viruses with implications for understanding recombination in SARS-CoV-2 evolution and the structure of clonal relationships and gene pools shaping its origins.

Core genes can have higher recombination rates than accessory genes within global microbial populations.

Recombination is essential to microbial evolution, and is involved in the spread of antibiotic resistance, antigenic variation, and adaptation to the host niche. However, assessing the impact of homologous recombination on accessory genes which are only present in a subset of strains of a given species remains challenging due to their complex phylogenetic relationships. Quantifying homologous recombination for accessory genes (which are important for niche-specific adaptations) in comparison to core genes (which are present in all strains and have essential functions) is critical to understanding how selection acts on variation to shape species diversity and genome structures of bacteria. Here, we apply a computationally efficient, non-phylogenetic approach to measure homologous recombination rates in the core and accessory genome using >100,000 whole genome sequences from Streptococcus pneumoniae and several additional species. By analyzing diverse sets of sequence clusters, we show that core genes often have higher recombination rates than accessory genes, and for some bacterial species the associated effect sizes for these differences are pronounced. In a subset of species, we find that gene frequency and homologous recombination rate are positively correlated. For S. pneumoniae and several additional species, we find that while the recombination rate is higher for the core genome, the mutational divergence is lower, indicating that divergence-based homologous recombination barriers could contribute to differences in recombination rates between the core and accessory genome. Homologous recombination may therefore play a key role in increasing the efficiency of selection in the most conserved parts of the genome.

Interplay of motility and polymer-driven depletion forces in the initial stages of bacterial aggregation.

Motile bacteria are often found in complex, polymer-rich environments in which microbes can aggregate via polymer-induced depletion forces. Bacterial aggregation has many biological implications; it can promote biofilm formation, upregulate virulence factors, and lead to quorum sensing. The steady state aggregation behavior of motile bacteria in polymer solutions has been well studied and shows that stronger depletion forces are required to aggregate motile bacteria as compared with their nonmotile analogs. However, no one has studied whether these same trends hold at the initial stages of aggregation. We use experiments and numerical calculations to investigate the polymer-induced depletion aggregation of motile Escherichia coli in polyethylene glycol solutions on short experimental timescales (∼10 min). Our work reveals that in the semi-dilute polymer concentration regime and at short timescales, in contrast to what is found at steady state, bacterial motility actually enhances aggregate formation by increasing the collision rate in viscous environments. These unexpected findings have implications for developing models of active matter, and for understanding bacterial aggregation in dynamic, biological environments, where the system may never reach steady state.

Food Polyelectrolytes Compress the Colonic Mucus Hydrogel by a Donnan Mechanism.

Systems consisting of a polyelectrolyte solution in contact with a cross-linked polyelectrolyte network are ubiquitous (e.g., biofilms, drug-delivering hydrogels, and mammalian extracellular matrices), yet the underlying physics governing these interactions is not well understood. Here, we find that carboxymethyl cellulose, a polyelectrolyte commonly found in processed foods and associated with inflammation and obesity, compresses the colonic mucus hydrogel (a key regulator of host-microbe interactions and a protective barrier) in mice. The extent of this polyelectrolyte-induced compression is enhanced by the degree of polymer negative charge. Through animal experiments and numerical calculations, we find that this phenomenon can be described by a Donnan mechanism. Further, the observed behavior can be quantitatively described by a simple, one-parameter model. This work suggests that polymer charge should be considered when developing food products because of its potential role in modulating the protective properties of colonic mucus.

High-molecular-weight polymers from dietary fiber drive aggregation of particulates in the murine small intestine.

The lumen of the small intestine (SI) is filled with particulates: microbes, therapeutic particles, and food granules. The structure of this particulate suspension could impact uptake of drugs and nutrients and the function of microorganisms; however, little is understood about how this suspension is re-structured as it transits the gut. Here, we demonstrate that particles spontaneously aggregate in SI luminal fluid ex vivo. We find that mucins and immunoglobulins are not required for aggregation. Instead, aggregation can be controlled using polymers from dietary fiber in a manner that is qualitatively consistent with polymer-induced depletion interactions, which do not require specific chemical interactions. Furthermore, we find that aggregation is tunable; by feeding mice dietary fibers of different molecular weights, we can control aggregation in SI luminal fluid. This work suggests that the molecular weight and concentration of dietary polymers play an underappreciated role in shaping the physicochemical environment of the gut. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Polymers in the gut compress the colonic mucus hydrogel.

Colonic mucus is a key biological hydrogel that protects the gut from infection and physical damage and mediates host-microbe interactions and drug delivery. However, little is known about how its structure is influenced by materials it comes into contact with regularly. For example, the gut abounds in polymers such as dietary fibers or administered therapeutics, yet whether such polymers interact with the mucus hydrogel, and if so, how, remains unclear. Although several biological processes have been identified as potential regulators of mucus structure, the polymeric composition of the gut environment has been ignored. Here, we demonstrate that gut polymers do in fact regulate mucus hydrogel structure, and that polymer-mucus interactions can be described using a thermodynamic model based on Flory-Huggins solution theory. We found that both dietary and therapeutic polymers dramatically compressed murine colonic mucus ex vivo and in vivo. This behavior depended strongly on both polymer concentration and molecular weight, in agreement with the predictions of our thermodynamic model. Moreover, exposure to polymer-rich luminal fluid from germ-free mice strongly compressed the mucus hydrogel, whereas exposure to luminal fluid from specific-pathogen-free mice-whose microbiota degrade gut polymers-did not; this suggests that gut microbes modulate mucus structure by degrading polymers. These findings highlight the role of mucus as a responsive biomaterial, and reveal a mechanism of mucus restructuring that must be integrated into the design and interpretation of studies involving therapeutic polymers, dietary fibers, and fiber-degrading gut microbes.

Target Turing patterns and growth dynamics in the chlorine dioxide-iodine-malonic acid reaction.

We study the growth dynamics of Turing patterns in the chlorine dioxide-iodine-malonic acid reaction-diffusion system in response to perturbations with visible light. We describe several mechanisms by which Turing patterns reappear after they are suppressed by illumination with a disc-shaped geometry. We observe that under specific conditions the patterns reorganize from a random configuration of spots and stripes to a set of ordered, concentric rings, which we refer to as target Turing patterns. These patterns closely resemble the unit cells of the Turing hexagonal superlattices known as black eye patterns. However, these target Turing patterns are not part of a larger superlattice structure, and they usually have a larger number of concentric rings. Numerical simulations support the experimental findings.